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A completely model-independent method of obtaining information on the spin using the quantum 
interference effect among various helicity states was proposed in a recent paper. Here we point out 
that this effect should be demonstrable in the existing data on e~e + — * W + W~ at LEP-II and 
pp — > Z° + j at Tevatron. 



There are many reasons to expect that new particle 
degrees of freedom will be discovered at the TeV en- 
ergy scale (Terascale), starting with the Large Hadron 
Collider (LHC) coming online later this year. The fact 
that the Terascale must have interesting physics has 
been known since Fermi's 1933 theory of nuclear beta 
decay which introduced a dimensionful constant Gf ~ 
(0.3 TeV) -1 . In its more modern incarnation, this con- 
stant represents the size of the Bose-Einstein condensate 
that makes the universe a gigantic superconductor. The 
analog of the Meissner effect then makes the range of the 
weak force as short as a billionth of a nanometer. 

At the least we expect the gap excitation of the super- 
conductor, the Higgs boson, to be discovered at the LHC. 
In addition, the quantum instability of this energy scale 
suggests new particles below a TeV in order to protect it 
from diverging to infinity. Many theoretical frameworks 
have been proposed in the literature: new strongly cou- 
pled gauge theory (technicolor [TJH]), fermionic dimen- 
sions of spacetime (supersymmetry [3]), bosonic dimen- 
sions of spacetime (extra dimensions [U [5] ) , new hidden 
extra symmetries (little Higgs [B]), Higgless theories [7J|8] 
etc. Many of these also provide candidates for the mys- 
terious dark matter of the universe. With great antici- 
pation the community awaits the imminent discovery of 
such exotic new particles in the upcoming LHC experi- 
ments. 

Once new particles are discovered, determining what 
theoretical framework they belong to is of foremost im- 
portance. For this purpose truly basic measurements will 
be required: mass, parity, and spin of the new particles. 
Among these, the spin measurement is both the key and 
the most challenging. Numerous studies exist that try to 
formulate strategics for spin measurements at the LHC 

[3|Tni[nm[T3m|TSl|TE]. Unfortunately, it is very dif- 
ficult to avoid model-dependent assumptions in the pro- 
posed measurement strategies. 

In a recent paper [T7], three of us (M.B., W.K., H.M.) 
proposed a completely model- independent way of obtain- 
ing information about spin at collider experiments. 1 The 
key clement is quantum interference among various helic- 



1 This possibility was originally suggested in |18l . 



ity states of the new particle, which, to our surprise, has 
not been discussed in the modern literature (see, how- 
ever, [HI]). We discussed how this method may work to 
discriminate the smuon in supersymmetry or the Kaluza- 
Klein muon in extra dimensions at the proposed Interna- 
tional Linear Collider (ILC). 

In this letter, we point out that the effectiveness of our 
proposed method should be demonstrable in the exist- 
ing data. In particular, e"e + — > W + W~ at LEP-II and 
pp — > Z a + j at Tevatron should allow highly significant 
studies of the quantum interference among helicities, and 
demonstrate the spin-one nature of the W and Z bosons 
without any model assumptions. As discussed in [17] . 
this method works particularly well close to the produc- 
tion threshold. This is good news for the LHC, as new 
physics there will likely be dominated by the energy range 
just above threshold. 

The proposed strategy is extremely simple. In order 
to obtain model-independent information about spin, or 
angular momentum in general, we resort to the general 
principles of quantum mechanics. The angular momen- 
tum operators generate spatial rotations; the unitary op- 
erator U(<j>) = e %J '^l n rotates space around the axis </> by 
the angle If we choose the rotation axis to be the mo- 
mentum vector of a free particle, it isolates the spin com- 
ponent because the orbital angular momentum is always 
orthogonal to the momentum vector L-p = (xxp)-p = 0. 
Therefore, the angular momentum along the momentum 
vector is nothing but its helicity, h — (s-p)/\p\. The rota- 
tion around the momentum axis by an angle <f> therefore 
gives the phase e lh ^ to the quantum mechanical ampli- 
tudes. 

Obviously a single phase factor does not lead to a phys- 
ical observable since the probability does not depend on 
phases. However, an interference effect may pick up the 
differences in phases among interfering amplitudes. For- 
tunately, particles produced in collisions are often in a 
linear superposition of various helicity states, which in- 
terfere when they decay into a common final state. This 
interference of different helicity states produces a cross 
section dependent on the coherent sum of individual ma- 
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trix elements squared: 



(TOC y] Alprod. (h)Md eca ,y(h, <j 
h 

M dccay (h, 0) = e lh *M dccay (h, <p = 0). 



(1) 



Here A4 pi odXh) an d M decay (h, = 0) are the production 
and decay matrix elements, which depend in detail on 
the helicity state h. However, all 4> dependence has been 
factored out into the exponential. It is clear from this 
sum that the azimuthal angular dependance of the event 
distributions N = a x £ (where C is the luminosity) is 



dN da 



x C = A + Ax cos ( 



A n cos(n4>), (2) 



where n — Ah is the difference between the highest and 
lowest helicity states contributing to the sum in Eq. ([I]). 
In this way, we obtain an unambiguous lower limit on the 
spin of the particle, s > (Ah)/2. As we will see, this limit 
is saturated, s = Ah/2, in the examples below, and the 
presence of the highest mode is clearly visible in collider 
data given sufficient statistics. 

In the cases of e~e + — » W + W~ with leptons plus jets 
final states and pp — > Z° + j with decays to electrons, 
spin-1 particles are produced in a superposition of helicity 
states. 2 In both cases, the event is fully reconstructable 
using the visible momentum in the event, and hence the 
angle 4> can be fully determined from data. 

The angle cf> is defined in the lab frame of the event as 
the angle between the production plane described by the 
W + W~ or Z°+j and the decay plane containing the lep- 
tonic decay products from the vector bosons. We define 
the positive z axis in the lab frame of LEP-II (Tevatron) 
as the direction of e~ (proton) beam, then the cosine of 
6 at LEP-II can be calculated as follows: 



^prod. 



Z X p w ± 



"decay 



p W ± X p e ± 



\z x p w ±\ " \Pw± x Pi±\ 

COS0 = "prod. ' "decay, (3) 



where p^± is the charged lepton from the decay of the 
W ± boson. The definition of <j> at Tevatron is the same 
as in Eq. ^ with the substitution of Z° for W ± . An 
arbitrary (but consistent) choice must be made to define 
which side of the production plane will contain positive 
4>. For LEP-II, we chose this positive direction to be 
in the direction of z crossed with the momentum of the 
lcptonically decaying W ± . Similarly, we chose the direc- 
tion of the proton beam crossed with that of the Z° at 
Tevatron (see Fig. [TJ. Based on our argument above, we 
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FIG. 1: The event kinematics of e e 

at LEP-II and pp -> Z° + j -> e"e + + j at Tevatron. The 
plane of pair produced vector bosons and the plane formed by 
the leptonic decay of one boson are shown. The angle <j> is the 
relative azimuthal angle between these two planes, defined in 
the lab frame of the event, as defined in Eq. <J3|) . Positive (f> 
are in the direction of the e~ (p) beam momentum crossed 
with the W~ (Z°) momentum for LEP-II (Tevatron). 



(GeV) 


c (pb- 1 ) 


182.25 


56.8±0.3 


188.63 


174.2L0.8 


191.58 


28.9±0.1 


195.52 


79.9±0.4 


199.52 


86.3±0.4 


201.62 


41.9±0.2 


204.86 


81.4L0.4 


206.53 


133.2L0.6 



TABLE I: LEP-II integrated luminosity £ as a function of 
beam energy yfs [20| . 



expect to see cross sections for these events as in Eq. ^ 
with n — 2. 3 

The LEP-II luminosity from the years 1997-2000 [20] 
are reported in Table [TJ The OPAL collaboration has ob- 
served 1574 events identified as qqev and an additional 
1573 qq[iv events [2JJ. Due to the low purity of the qqrv 
sample, we ignore those events. Similar data sets are 
available from the ALEPH [20J, DELPHI [22], and L3 
[25] collaborations. 

The CDF collaboration has data for Z° + j consisting 
of 6203 events [21] after selection cuts from 1.7 fb _1 of 
luminosity at 1.96 TeV beam energy. D0 has a similar 
data set available [25]. A total luminosity of 8 fb _1 is 
expected to be available from Tevatron at the conclusion 
of data collection. 

Parton level matrix elements for W + W~ and Z° + j 
(where the jet consists of a gluon or first generation (anti) 
quark at the parton level) production were calculated 



2 It is for this reason we cannot consider pp — > Z without jets. 
In such events, the Z is produced in only one spin state, de- 
pending on the spin of the initial state quarks. While the cross 
section would contain a sum over Z helicity, the sum would be 
incoherent. 



3 It should be noted that if the collider beams are identical, this 
choice of positive <j> suffers from an ambiguity which maps <f> — > 
<p + n. This may, for example, introduce difficulties in measuring 
A n (n odd) parameters at LHC. 
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in HELAS [25] , while the numerical integration program 
BASES [27] was used to determine the differential cross 
section and integrate over all other kinematic variables. 
For the simulation of the Tevatron results, a K factor 
of 1.4 was used to correct for higher order QCD effects, 
in accordance with [21], and CTEQ5L PDFs were imple- 
mented using LHAPDF [28 . The Tevatron results and 
fits were confirmed using ALPGEN [29 . 
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FIG. 2: Differential distribution of events do/d(f> x C for a) 
e~e + —> W + W~ — > qqt^v using the bEP-II run data in 
Table [i] and b) pp — > Z° + j — > e e + + j with luminosity 
C = 1.7 fb _1 . No cuts are applied on the LEP-II simulation, 
Tevatron results have pt > 30 GeV and |r/| < 2.1 on the jet. 





LEP-II 


Tevatron 


Ai/A 


-0.267L0.023 


0.036±0.009 


A 2 /A 


-0.085±0.025 


0.100±0.009 


A s /A 


0.000±0.025 


0.000±0.009 


A 4 /A 


0.000±0.026 


0.000±0.010 



TABLE II: Fits to the parameters A n in Eq. (Pif for the dif- 
ferential distributions of e"e+ -> W + W~ -> qqi^v (LEP-II) 
using the integrated luminosity in Table]!] and pp — > Z° +j — > 
e~e + + j (Tevatron) using C = 1.7 fb~ . Errors for each pa- 
rameter are obtained by marginalizing over the other four 
parameters in the fit. No cuts are applied on the LEP-II sim- 
ulation, Tevatron results have pt > 30 GeV and I77I < 2.1 on 
the jet. 

The generated histograms are assigned Gaussian sta- 
tistical error bars based on the realistic experimental lu- 
minosities. However, no statistical fluctuations are as- 
signed to the central values. As a consequence, the fit 
results correspond to an average experiment |30j . 

Before the application of cuts, the differential cross 
sections for the two processes of interest show a clear 
cos 4> and cos 2<j> dependence, as expected for the de- 
cays of spin-1 bosons. These distributions are shown in 
Fig. [2] We then fit the parameters Aq,A\,A2,A^, and 
A4 in Eq. ^ to the event distributions. 4 For each of the 
five parameters A n , 1-er error bars are calculated after 
marginalizing over the other four. Results for the LEP- 
II and Tevatron simulations are shown in Table In] in 



Jet transverse momentum 


Pt,j > 30 GeV 


Jet r\ 


hi < 2.1 


Invariant mass of lepton pair 


66 < mu < 116 GeV 


Central electron 77 


77 < 1 


Second electron 77 


\rj\ < 1 or 1.2 < |r?| < 2.8 


Electron Et 


E T > 25 GeV 


Electron isolation cuts 


AR ej > 0.7 



TABLE III: Event selection cuts imposed by the CDF collab- 
oration on pp — * Z° + j — > e~e + + j events. In each event, 
one electron must be central, and pass stricter cuts than the 
second electron. The isolation cut parameter is defined as 
y/ (A</>) 2 + (A77) 2 = AR [H]. 



Lepton momentum 


pi > 25 GeV 


Polar angle 9 of final state particles 


| cos 6»| < 0.95 


Neutrino energy fraction 


R v > 0.07 


Visible energy fraction 


i? vis > 0.3 


Neutrino transverse momentum 


Pt,„ > 16 GeV 


Lepton isolation 


AR > 0.75,0.5,0.2 



4 These fits are to the numerically integrated differential cross- 
section, not generated events. 



TABLE IV: Event selection cuts imposed by the OPAL col- 
laboration on e~e + — > W + W~ — > qql^u events. Energy 
fraction is defined as R a = E a /yfs, where a is either the neu- 
trino v or the total visible energy. The lepton isolation cut 
was implemented using x/(A<j)) 2 + (Arj) 2 = AR with a range 
of AR values rather than limiting total energy deposited in 
cone surrounding the lepton as in [3 1J . 



order to compare simulations with different numbers of 
events, values of A n /A are reported rather than A n . It 
it clear at this stage that the results are consistent with 
the decay of spin-1 bosons. 

However, cuts must be applied to the events recorded 
at LEP-II and Tevatron, both due to detector geometry 
and in order to reduce background. These cuts will affect 
the azimuthal distribution present in da /d<j) x C, and so 
can obscure the signal necessary for spin measurements. 
The Tevatron cuts (Table III) were taken from the CDF 
experiment [24], while the OPAL [31] cuts (Table [TVl 
were used to simulate the LEP-II data. 

Our simulation did not include parton showers or 
hadronization, so we could not implement the lepton iso- 
lation cut used by OPAL, which placed a limit on the to- 
tal energy deposited in a cone centered on the lepton. In- 
stead, we used a cut on AR = \J (Acf)) 2 + {A-q) 2 between 
the jet and the leptons. Three values for AR were used: 
0.2, 0.5 and 0.75, which gave total efficiencies for the 
cuts of 79%, 76%, and 72% respectively. The cuts used 
by the OPAL collaboration had an efficiency of 85% for 
final states with an electron and 89% for muons. The dis- 
tributions of the Tevatron and LEP-II (with AR = 0.75) 
simulations after cuts are shown in Fig. [4] 

Fitting the distributions to Eq. Q , we find the results 
in Table [V] These results clearly show that the imposed 
cuts introduce spurious high frequency modes. The cor- 
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LEP-II 




AR = 0.75 


AR = 0.5 


AR = 0.2 


Ax/ An 


-0.082±0.025 


-0.082±0.026 


-0.082±0.025 


A 2 /A 


-0.293±0.026 


-0.302±0.027 


-0.308±0.026 


A3/A0 


0.110±0.027 


0.114±0.028 


0.117±0.028 


A4/A0 


-0.099±0.028 


-0.099±0.029 


-0.096±0.029 






Tevatron 




Ai/Ao 


0.029±0.012 


A 2 /A 


-0.277±0.012 


A 3 /A 


-0.021±0.013 


A 4 /A 


-0.123±0.014 



TABLE V: Fits of the differential distribution of e"e+ -> 
W + W -> qq£ ± v (LEP-II) with the cuts in Table |iv] and 
pp -> Z°+j -> ri + +j (Tevatron) with the cuts in Table |ih| 
to parameters A n in Eq. pj. Luminosities are as in Table]!!] 
1-a errors for each parameter are obtained by marginalizing 
over the other four parameters in the fit. 



responding non-zero A 3 and A4 components may naively 
be confused for evidence of spin-2 particles. However, the 
cuts are responsible for introducing new <f> dependence by 
selecting out new directions relative to the production 
axis of the gauge bosons. 

We illustrate this effect for the case of cuts in the for- 
ward direction (large \rj\ and |cos0|) in Fig. [3] Here we 
see two decays which are kinematically identical in the 
boson rest frame save for azimuthal rotations. In Fig. [3^, 
the event survives the cuts, as neither lepton lies suffi- 
ciently close to the z axis. However, in Fig. [3J3, rotating 
the decay plane about the axis of the boson momentum 
yields an event which is eliminated by the cuts. This is 
the source of unwanted <f> dependences in the differen- 
tial distributions with cuts. Similar problems arise due 
to isolation cuts, which depend on the proximity of the 
leptons to the other particles in the final state, as well as 
cuts on leptonic transverse momentum. 

Since this <f> dependence did not arise from the quan- 
tum interference of helicity amplitudes, we cannot expect 
the <fi dependence of the cross section to accurately reflect 
the spin of the decaying particles. Thus non-zero A3 and 
A4 components do not indicate a higher spin state, but 
rather a breakdown of the proposed spin-measurement 
technique. 

The solution to this problem is relatively straightfor- 
ward. For new azimuthal dependences to be avoided, the 
cuts cannot pick out 'special' directions relative to the 
original momentum of the decaying boson. Therefore we 
impose 'rotationally invariant cuts' in which we require 
that each event not only passes the experimental cuts 
but continues to do so when the decay plane is rotated 
around the boson production axis. This avoids the intro- 
duction of a new directional dependence since we restrict 
ourselves to only those events which could never overlap 
the forbidden regions of the detector regardless of orien- 
tation. However, these cuts are very inefficient: the cuts 
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FIG. 3: A depiction of the detector volume demonstrating the 
rotational dependence induced by the cuts. The shaded for- 
ward regions (large values of \r)\ and | cos#|) are inaccessible 
due to detector geometry and background cuts. Two sample 
events are depicted in a) and b). These events are kinemat- 
ically identical in the boson rest frame save for a rotation in 
4>. The event a) survives the cuts, while the event b) fails. 
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FIG. 4: Differential distributions for a) e~e + — > W + W~ — > 

and A_R = 0.75 and b) pp — > 
from Table III Luminosities 



qqi v with the cuts in Table 



IV 



Z° +j^e e + +j with the cuts 
are as in Fig. [2] 



on LEP-II data preserve only 12% of the original events, 
while the cuts for the Tevatron leave less than 1%. 

The CDF cuts are very inefficient due to the small al- 
lowed 1 77 1 region for the central electron (see Table III I. 
Recent preliminary CDF measurements have demon- 
strated that the cuts can be relaxed while still main- 
taining a background level of less than 5% [31] ■ These 



loosened cuts are identical to those in Table III for pt 
and 77 of the jet and the invariant mass of mu- However, 
the central lepton is allowed E T > 20 GeV and \t)\ < 2.6, 
while the second electron must have Et > 10 GeV and 
\ri\ < 2.6. If both leptons have 2.6 > \r]\ > 1.0, E T must 
be greater than 25 GeV. Finally, AR e j must be greater 
than 0.4. With these relaxed numbers, the total num- 
ber of events in the simulated sample is 5821 and the 
efficiency of the rotationally invariant cuts is 18%. 
The result of these rotationally invariant cuts on the 
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FIG. 5: Differential distributions for a) e~e + — > W + W~ — > 
qql^v and b) pp — > ^° +j — > e~e + +j requiring rotationally 
invariant cuts. Luminosities are as in Fig. [2] 



LEP-II and Tevatron data are shown in Fig. [5] (compare 
to Fig. [2| . Table VI confirms that this technique restores 
the 4> dependence expected by the interference argument. 



In the case of the Tevatron results with loosened cuts, 
the data is clearly consistent with the Z being a spin-1 
vector boson. The A\ parameter is non-zero at 1.8a, the 
A2 parameter is non-zero at nearly 4<r, and the higher 
modes are consistent with zero. It is important to re- 
call that a lower bound on the spin is obtained from the 
highest non-zero mode, therefore the Act signal in Ai is 
far more important then the 1.8cr deviation from zero in 





LEP-II 




Ai? = 0.75 


AR = 0.5 


Ai? = 0.2 


A1/A0 


-0.215±0.069 


-0.214±0.060 


-0.207T0.053 


A2/A0 


-0.068±0.071 


-0.071±0.062 


-0.072T0.055 


A 3 /A 


0.000±0.073 


0.000±0.064 


0.000±0.057 


A 4 /A 


0.000±0.075 


0.000±0.065 


0.000±0.058 






Tevatron 




Ai/Ao 


0.039±0.022 


A 2 /A 


0.083±0.021 


A 3 /A 


0.000±0.022 


Ai/Ao 


0.000±0.023 



TABLE VI: Fits of the differential distribution of e~e + -> 
W+W -> qql^v (LEP-II) and pp -> Z° + j -> t~l+ + j 
(Tevatron) to the parameters A„ in Eq. pj, requiring events 
that pass the cuts in Tables IV and |III| (with relaxed Er, \rj\, 
and Ai? cuts as described in the text) after rotation about 
the momentum axis of the decaying vector boson. The lumi- 
nosities are the same as in Tables HIl and IVl 1-a errors for each 
parameter are obtained by marginalizing over the other four 
parameters in the fit. 



and increase our confidence in the result correspondingly. 
Using, for example, the estimated total integrated lumi- 
nosity of 8 flV 1 for the Tevatron, the parameters have 
the values shown in Table |VII| Another possibility is to 
use the muon decays of the Z u . However, the rotationally 
invariant cuts will likely take a high toll on such events, 
as the muon tracking system at CDF extends only up to 

h = 1.5 im- 

The situation with the LEP-II simulation is more com- 
plicated. While the Ai parameters are non-zero at over 
3cr, the A2 parameters differ from zero by only one stan- 
dard deviation. A larger data set would of course solve 
this problem. As all four LEP-II experiments (ALEPH, 
DELPHI, L3, and OPAL) have approximately equal 
statistics available, a two-fold increase in the statistical 
significance could be achieved by combining the events 
from these collaborations; the resulting ratios A n j Aq are 
shown in Table IVTT1 

Another possibility is that some reduction in required 
cuts would increase the efficiency of the rotationally in- 
variant cuts without greatly degrading the sample pu- 
rity. A likely candidate for this in our analysis is the Ai? 
cut, which was introduced as a stop-gap measure to ap- 
proximate the jct-lcpton proximity cut used in the OPAL 
analysis. However, even with the value of A R = 0.2, the 
efficiency of the cut is lower than the 85% reported by 
OPAL. Setting Ai? = is clearly an unrealistic cut, but 
as demonstrated in Table |VII| indicates the possibilities 
offered by higher statistics. 

In conclusion, we have demonstrated that the quan- 
tum interference among the matrix elements of different 
helicity states provides model-independent probe of par- 
ticle spin. Using realistic data sets, rotationally invariant 
cuts can be implemented which correct for the spurious 
high-frequency noise introduced by the cuts imposed by 
detector geometry and background reduction. Though 
these techniques come at a price in terms of efficiency, it 
seems possible to relax the cuts in such a way that the 
weak gauge boson spins can be measured at sufficient 
significance at current colliders. 

Measurements of the spin of new particles is expected 
to be a critical discriminator of new physics at the LHC. 
As a result, techniques such as the one proposed here 
are very important. Though the spins of the W and Z 
bosons are not in doubt, we find it encouraging that this 
new method can be tested on the available data. Such 
work would be of great use in the coming LHC era. 



From these results there is always the possibility that 
the parent Z is a higher spin particle and that some con- 
spiracy amongst the matrix elements in Eq. (|T|) prevents 
the A3 and A4 terms from appearing in the sum. In this 
interpretation, we can still state unambiguously that the 
Z is at least spin-1, and that the data suggest it is not of 
higher spin. 

Higher statistics would allow a reduction of error bars 
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LEP-II 


Tevatron 




Combined 


AR = 


C = 8 fb" 1 


Ai/A 


-0.207±0.027 


-0.211±0.050 


0.039±0.010 


A 2 /Ao 


-0.072±0.028 


-0.081±0.052 


0.083±0.010 


A 3 /A 


0.000±0.028 


0.000±0.053 


0.000±0.010 


A4/A0 


0.000±0.029 


0.000±0.054 


0.000±0.010 



TABLE VII: Fits of the differential distribution to the pa- 
rameters An in Eq. (2) for: e~e+ -> W + W~ -> qqt±v with 
the jet-lepton cut parameter A_R = 0.2 and combining the 
data sets of ALEPH, DELPHI, L3, and OPAL (LEP-II, Com- 
bined), e"e+ -> W + W~ -> qql^v with the OPAL data set 
and AR set to zero (LEP-II, AR = 0), and pp -> Z° + j - -> 
£~£ + +j with 8 fb _1 integrated luminosity (Tevatron). We 
require all ev ents pass the cuts in Tables IV (with AR as in- 
dicated) and III (with relaxed Et and \rj\ cuts as described 
in the text) after rotation about the momentum axis of the 
decaying vector boson. 1-a errors for each parameter are ob- 
tained by marginalizing over the other four parameters in the 
fit. 
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